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ABSTRACT 

The inclination of M3 1 is too close to edge-on for a bar component to be easily recognised 
and is not sufficiently edge-on for a boxy /peanut bulge to protrude clearly out of the equatorial 
plane. Nevertheless, a sufficient number of clues allow us to argue that this galaxy is barred. 
We use fully self-consistent iV-body simulations of barred galaxies and compare them with 
both photometric and kinematic observational data for M31. In particular, we rely on the 
near infrared photometry presented in a companion paper. We compare isodensity contours 
to isophotal contours and the light profile along cuts parallel to the galaxy major axis and 
offset towards the North, or the South, to mass profiles along similar cuts on the model. All 
these comparisons, as well as position velocity diagrams for the gaseous component, give 
us strong arguments that M31 is barred. We compare four fiducial TV -body models to the 
data and thus set constraints on the parameters of the M31 bar, as its strength, length and 
orientation. Our 'best' models, although not meant to be exact models of M31, reproduce in 
a very satisfactory way the main relevant observations. We present arguments that M31 has 
both a classical and a boxy/peanut bulge. Its pseudo-ring-like structure at roughly 50' is near 
the outer Lindblad resonance of the bar and could thus be an outer ring, as often observed 
in barred galaxies. The shape of the isophotes also argues that the vertically thin part of the 
M31 bar extends considerably further out than its boxy bulge, i.e. that the boxy bulge is only 
part of the bar, thus confirming predictions from orbital structure studies and from previous 
iV-body simulations. It seems very likely that the backbone of M3 1 's boxy bulge is families of 
periodic orbits, members of the Xi-tree and bifurcating from the Xi family at its higher order 
vertical resonances, such as the xlv3 or xlv4 families. 



1 INTRODUCTION 



It is very easy to recognise bars in disc galaxies which 
are far from edge-on. This, however, is not true for discs 
seen edge-on, the safest appro a ch being to use k i nemat- 
ics iKuiiken & Memfield Il995t iMerrifield & Kuiikenl | l99S : 
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Athanassoula &BureaJll999T Chuns & Bureau 2004). Yet kine- 
matics are often not available, and thus one has to rely on mor- 
phology and photome try alone. A large number of properties (see 
lAthanassoulal l20"05al) . hereafter A05, for a global discussion and a 
compilation of references) argue that a boxy/peanut bulge is a part 
of a bar seen edge-on. Thus, the mere presence of su ch a bulge is 
telltale of the existence of a bar component (but see IPatsis et al. I 
120021) . arguing that other rotating bisymmetric perturbations could 
also form boxes). If the galaxy is observed side-on (i.e. edge-on 
with the line of sight along the bar minor axis), then the boxy or 
peanut feature is clearly visible. If, however, the bar is seen end- 
on (i.e. edge-on with the line of sight along the bar major axis), 
then the bar has the outline of a classical bulge and it can not 

1 Following A05, we distinguish throughout this paper three types of 
bulges. Classical bulges were formed in the early stages of the galaxy for- 
mation process, before the present discs, from gravitational collapse, or hi- 
erarchical merging of smaller objects and the corresponding dissipative pro- 



be recognised. Luckily, even deviatio ns as small as 10° from end- 
on w ill reveal the bar/peanut feature iLutticke. Dettmar & Pohler] 
l200d A05). Evidence for bars in edge-on sys tems can also be found 
from their photometry iLiitticke et al. l200u) . Indeed, in such cases, 
cuts along the equatorial plane reveal plateaux on either side of the 
center. Such plateaux are also seen in similar cuts in bar red galaxy 
simulations (Fig. 6 in lAthanassoula & Misiriotisl 120021) . hereafter 
AM02) where it can be verified that they extend roughly to the 
end of the bar. The existence of such plateaux on equatorial cuts 
is, however, not unambiguous evidence for the existence of a bar 
component, since these features can also be due to lenses. 

Although considerable effort has been put into recognising 
bars in e dge-on systems, very littl e has b een done for near-edge-on 
systems. iBureau & Athanassoulal 120051) have described the posi- 
tion velocity diagrams (hereafter PVDs) of simulations viewed un- 
der such an orientation, as well as the corresponding V, a, hs and 
hi profiles (see their Fig. 4). If, however, kinematical data are not 



cesses. Boxy/peanut bulges are part of the bar, i.e. they are constituted by 
stellar disc material, and they form naturally from the long-term evolution 
of the stellar component of bar unstable discs, as witnessed in a very large 
number of simulations. Finally, disc-like, or discy bulges form out of gas 
that concentrates to the inner parts of the disc under the influence of the 
gravitational torque of a bar, and subsequently forms stars. 
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available, then recognising a bar in such orientations is more diffi- 
cult than in edge-on galaxies. Indeed, in such cases the boxy/peanut 
bulge does not clearly stick out of the equatorial plane, so that this 
diagnostic can not be used. Furthermore, no photometric analysis 
of simulations viewed under such orientations is available to pre- 
dict what one should expect for the surface density profiles. 

Although many disc galaxies are viewed at angles near to 
edge-on and a large fraction of them should have boxy/peanut 
bulges and bars, very few have been discussed in the literature. 
NGC 7582 has an inclinat ion of approximatel y 65° and, when 
viewed in the near infrared JOuillen et fl/.ll997l) . allows for visual 
identification of both the bar and the peanut. Simila r comments can 
be made for NGC 4442 iBettoni & Gallettall 19941 . viewed at ap- 
proximately 72°. M31, our nearest big neighbour, is near to edge- 
on and exhibits isophotal shapes similar to those of NGC 7582 and 
NGC 4442. Could it be barred? We will aim to answer this ques- 
tion by comparing results from iV-body simulations to the relevant 
data. 

M31 (NGC 224) is an Sb type spiral galaxy, which has an in- 
clination angle of about 77°. The isophotes of M3 1 have a twist 
in the inner parts jLindblaj|lj56yHoj ge&Kennicutt|[l^ and 
in the outer parts (Walterbos & Kennicutll988t) . the latter presum- 
ably due to a warp. In the main region of the disc, however, the 
position angle is well defi ned and nearly constant (see Fig. 7 of 
iHodee & Kennicuttlll982l) . We will thus henceforth ad opt a fidu- 
cial disc position angle of 38° IdeVaucouleursI ll958f) and later 
work). Further in, in the region often referred t o as the 'bulge re- 
gion' , the position angle is 10° to 20° larger faichter & Hognej 
( 1963) a nd later work). The distance of M31 is estimated to be 
783 kpc lHollandll998l) . which makes 10' roughly equal to 2.28 
kpc. In earlier studies, however, it was assumed that 10' is roughly 
equivalent to 2 kpc. For this reason, and in order to avoid confu- 
sion, we will here use arc minutes to measure distances. Due to 
its proximity and its interesting structure, M31 has been the sub- 
ject of a very large number of studies, of which a complete list of 
references is beyond the scope of this paper. A mongst the older 
photome try, we have used isop hotal results from lLindblad ll956t) 
and fromlde y aucouleursNl958l) , as well as position angle measure- 
ments froin lHod^e^^Kermicutl ll982l) . More recently, M31 was 
observed by the 2MASS facility as part of the "6X Survey", going 1 
magnitude deeper than the original 2MASS survey 2 . These images, 
in J, H and K 3 were analysed by Beaton et al. (2006, hereafter Pa- 
per II) revealing important morphological properties, especially of 
the bulge. It is these observational results that originally motivated 
the work described here. 

The id ea that M31 is b arred is not new. It was initially pro- 
posed bv B. lLindbladNl956l) . based on the analysis of the isophotal 
shapes of a deep red exposure. Approximating the isophotes by el- 
lipses, he found that, on the plane of the sky, their major axis is 
at an angle of roughly 10° south of the galaxy major axis and that 
the average projected axial ratio is 1.6. He proposed that the bulge 
is in fact a triaxial bar with axial ratios 1:0.56:0.28 or 1:0.56:0.36 
and semi-major axis length around 15'. Yet, the existence of a bar 
in M31 was generally not accepted and soon forgotten, so that 
hardly any of the ve ry numerous subsequent papers on M3 1 men- 
tion this possibility. IStarta ll977h returned to this conclusion with 
his calculation of a one-parameter family of triaxial models for the 
M31 bulge. Their axial ratio in the galaxy plane is about 1.6 to 



2 See the 6X description at 

http://www.ipac.caltech.edu/2mass/releases/allsky/doc/ 



1.8, depending on the angle of the bar major axis w ith the gal axy 
major axis. The two solutions initially proposed by iLindblad are 
mem bers of thi s famil y. To distinguish between IStarkh models, 
IStark & Binnevl (1994) modeled the gas flow in the inner couple 
of arc minutes using a Ferrers ellipsoid and a spherical component. 
They propose that, measured on the plane of the galaxy, the angle 
between the bar semi-major axis and the projection of the line of 
sight onto the plane of M3 1 is about 20° . 

In this paper we will examine the possibility of M31 being 
a barred galaxy based partly on the photomet ric results of pa- 
per II and partly on observed kinematics {brinks & Shandfl984l 
iBrinks & Burtor]ll984l : lRubin & Fordll970n . The wealth of infor- 
mation provided by these results will allow us to go further than 
just debate on the possible existence of a bar. We will thus aim to 
find some of the main properties of this bar. It should, however, 
be clear that we will not present a complete model of M31, which 
would be well beyond the scope of the present paper, but we will 
instead try to investigate all the available observational clues in or- 
der to characterize and constrain the bar properties. This paper is 
organised as follows. Section|2|describes the simulations and Sec- 
tion[3]compares them to the photometrical observations of Paper II. 
The kinematics are presented in Section|4] We discuss our results 
in Section|5|and briefly summarise in Section|6| 



2 SIMULATIONS 

Both observed and simulated bars come in a large variety of 
lengths, strengths and shapes. We will thus not restrain ourselves 
to one single simulation to compare with the observations, but use 
four distinct cases, hoping to be able to set constraints on the pa- 
rameters of the bar 3 . 

AM02 introduced and described two different types of bars, 
which they termed MD and MH, respectively. MD-type bars grow 
in galaxies in which the disc dominates the dynamics in the inner 
parts, i.e. are maximum disc galaxies with haloes with extended 
cores. On the other hand, MH-type bars grow in galaxies in which 
the halo has a small core, i.e. is centrally concentrated, so that the 
halo contributes as much as, if not more than, the disc to the dynam- 
ics in the bar region. The bars growing in these two types of galax- 
ies have very different morphology and kinematics (AM02). MH- 
type bars are thinner and longer than MD types and their outline 
is rectangular-like, contrary to MD bars which have an elliptical- 
like outline. The amplitude of their m = 4, 6 and even 8 Fourier 
components of the density reaches a considerable fraction of the 
amplitude of their m = 2, contrary to MD-types in which the m = 
6 and 8 are negligible. The density profile along the bar major axis 
(face-on view) also differs in the two types of bars. In MH-types it 
has two plateaux, one on each side of the center, with abrupt drops 
at the end of the bar, while in MD-types it drops near-exponentially 
with distance from the center (Fig. 5 in AM02). Bars in MH-type 
models often have ansae and/or an inner ring, which is elongated, 
but not far from circular and has the sa me major ax is as the bar, as 
inner rings observed in barred galaxies <Butil986l) . Their side-on 

3 All simul ations described here have live haloes. This is necessary since, 
as shown by Athanassoula 1 2002^), in order for bars to evolve, the halo must 
be able to participate in the angular momentum exchange process and to 
respond to the disc evolution. Thus, in the simulations the halo must be 
live, i.e. composed of particles. Rigid haloes, i.e. haloes which are simply 
a non-evolving potential imposed on the disc particles, artificially hinder or 
quench bar formation. 
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shape evolves first to boxy and then to a peanut or 'X' shape, in 
contrast to the MD-types for which the side-on outline stays boxy. 
The side-on velocity field of MH-types shows cylindrical rotation, 
while that of the MD-types does not. More information on these 
properties can be found in AM02. There is a clear connection be- 
tween the bar morphology and the amount of angular momentum 
emitted by near-resonant material in the bar region and absorbed 
by near-resonant material in the halo, in the sense that more angu- 
lar mo mentum has been exch anged in MH-types than in MD-type 
models \ Athanas soulal2005bT) . Furthermore, correlations have been 
found between the bar str ength and the amoun t of angular momen- 
tum absorbed by the halo lAthanassoulefeOolh . 

Since MD and MH-type bars are so distinct, we decided to 
include one of each in our fiducial models. Our other two fiducial 
bars are also one MH and one MD type, but with a classical bulge 
component as well. In the notation of AM02, they are MDB and 
MHB types. It should be stressed at this point that, both in observa- 
tions and in simulations, there is a continuum of bar types ranging 
between those described above. Furthermore, other parameters than 
the halo mass distribution influence the amount of angular momen- 
tum exchanged within a galaxy and therefore its bar strength. These 
include the responsiveness of the halo (in part icular its veloc ity dis- 
persion), the velocity dispersion in the disc lToomre]| 19641 the O 
param eter) and the existence of a bulge component lAthanassoulal 
12003) . Thus, it is preferable to consider our four fiducial models as 
examples of bar types, rather than as absolute measures of the halo 
mass distribution. For simplicity and in order to facilitate links with 
previous work, we will refer to our four fiducial models as MD, 
MH, MDB and MHB. 

The initi al conditions of our four models are as described in 
AM02 and in lAthanassoulal 1200 3l) . Thus, the radial density distri- 
bution of the disc is exponential with mass Md and scalelength h 
and its vertical distribution follows a seek 2 law with scaleheight 
zq. Its radial veloc ity dispersion is given by the Toomre Q parame- 
ter | Tpomrell964) . The halo is described by eq. (2.2) of lHernauistl 
1 1993) and has two characteristic radii, 7 and r c , out of which the 
first one measures the core size and the second one the outer trun- 
cation . The bulge is described by a Hernquist sphere iHernauisl 
Il99d) of mass Mb and scalelength a. All the models described 
here have Md = 1, h = 1, Zo = 0.2, Mh = 5 and r c = 10. For 
the other parameters we have Q = (1.2, 1.2, 1, 1), 7 = (5, 0.5, 5, 
0.5), M b = (0, 0, 0.6, 0.4) and a = (-, -, 0.6, 0.4), for models (MD, 
MH, MDB, MHB), respectively. Simulat ions were run using W. 
Dehn en's treecode |DehnerJl200d. l2002al) . as described in AM02 
and in lAthanassoulal 12*003^ 

In these simulation there are about 10 6 particles in the halo 
component and between 200 000 and 400 000 particles in the vis- 
ible components (disc and classical bulge). These are amply suffi- 
cient for following the evolution of the run and the formation of the 
bar and boxy bulge. For accurate comparisons with the observed 
isophotes, however, one needs isodensities which have less noise 
than what is obtained with this number of particles. Rather than 
use some smoothing, which could also smooth out real features, we 
used the technique described in A05, which allows an a posteriori 
increase the number of particles. Namely, we considered ten snap- 
shots, closely spaced in time so that the bar will not have evolved 
much in between any two, except of course for its figure rotation. 
The shapshots were rotated, so that the major axis of the bar is at 
the same position angle in all ten cases, and then stacked. This in- 
creases the number of particles tenfold. We then assumed two-fold 
symmetry, i.e. symmetry with rotation by 180° and with respect to 
the equatorial plane. This allowed us to avoid asymmetries incon- 



sistent with M31 4 and which would, therefore, make all compar- 
isons more difficult. Thus, particle numbers were increased by a 
factor of 40 at the expense of a considerable amount of work be- 
fore each display. This technique allows us to achieve isodensity 
contours with very little noise without applying any convolution 
which could have smoothed out characteristics of the isodensities 
as well as the noise. 

Face-on and side-on views of the luminous material (i.e. the 
disc and, for MHB and MDB, also the classical bulge) can be seen, 
respectively, in the lower and the upper panels of Fig. Q respec- 
tively. Comparison of the face-on views shows that the MD bar 
is the fattest and shortest of the four and that its isodensities are 
elliptical-like. Model MH has the longest and strongest bar and its 
isodensities are rectangular-like. For the other two, the central parts 
of the isodensities, where the classical bulge contribution is very 
strong, are elliptical-like and the parts further out, where the bar 
dominates, are rectangular-like. This is simila r to what is observed 
in ma ny early type strongly barred galaxies lAthanassoula et al\ 
Il99d) . The MH bar is surrounded by an inner ring and the MHB 
bar has clear ansae. 

The side-on views (upper panels) also present many differ- 
ences. Model MD has a boxy-like shape, while MH displays a 
strong peanut, or 'X'-like feature. It was shown in A05 that in cases 
with bulges the classical bulge material 'fills' up the central and off 
the equatorial plane parts, so that the outline becomes boxy-like. 
Indeed, models MHB and MDB have boxy-like shapes. Structures 
as those seen in these models are referred to as box y/peanut b ulges 
and h ave been observed in many edge-on galaxies iLiitticke et al. I 
l200d) and their link to bars is by now well established (A05 and 
references therein). 

Since we can view our simulations both face-on and side- 
on, it is possible to obtain estimates of both the bar and of the 
boxy/peanut bulge extent. As discussed in AM02, there are several 
methods for doing this and their results show often considerable 
differences. For this reason, we adopted here very simple eye esti- 
mates. We obtained the bar length from the face-on view by finding 
by eye the extent of the last isophote which outlines the bar. The 
boxy bulge length was measured as a distance from the center and 
along the major axis on the side-on view. Within the boxy bulge ex- 
tent, the distance of the isodensities from the equatorial plane either 
stays roughly constant (boxy bulges), or increases with increasing 
distance from the center (peanut, or 'X' shape bulges). Beyond the 
boxy bulge extent, the distance of the isodensities from the equato- 
rial plane decreases with increasing distance from the center. Obvi- 
ously there is some uncertainty in these determinations. However, 
if sufficient care is taken, these measurements can be as reliable as 
those based on other methods. The bar and the boxy/peanut bulge 
extent are given in Fig. Q by the solid and dashed lines, respec- 
tively, and have been extended through all three panels to allow us 
to compare the two extents. Our figure shows clearly that the size 
of this boxy bulge is shorter than the bar length and that the part of 
the bar outside the boxy bulge region is vertically thin. This is in 
agree ment with previous simulations and w ith orbital structure re- 
sults IPatsis. Skokos & Athanassoulal2002T) and has been discussed 
extensively in A05. 

For model MH we plotted two solid lines, i.e. we give two 
possible estimates of the bar length. Indeed, if we consider that the 
bar stops where the ring starts, then the shortest of the two bar- 



4 The asymmetries in M3 1 are presumably linked to its two main compan- 
ions, M32 and NGC 205, which are not included in the modeling here. 
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Figure 1. Three views of the simulations which will be compared to M3 1 . Simulations MD, MH, MDB and MHB are shown in columns 1 to 4, respectively. 
The upper panels give the edge-on views and the lower one the face-on views. The middle ones are views at 77°, the inclination angle of the M31 disc. The 
grayscale is based on the logarithm of the projected density. The isodensities are chosen at levels that show best the shapes discussed in section lTTI The solid 
vertical line(s) in each row give(s) a rough estimate of the location of the end of the bar, as obtained from the face-on view. The dashed lines give a rough 
estimate of the extent of the vertically extended boxy feature (i.e. the length of the boxy bulge) as seen from the edge-on view. Note in the middle panel that 
isophotes whose major axis is within the dashed line are boxy-like. On the other hand, isophotes between the dashed and the full vertical lines have a pointed 
feature, elongated roughly in the direction of the bar. This is the signature of the vertically thin part of the bar. 



lengths is the appropriate one. Alternatively, if we consider that the 
bar continues until the end of the ring, then it is the longest of the 
two barlengths that is appropriate. The two alternatives will be dis- 
cussed further in Section l3~2l 

The middle row in Fig. \\\ gives an intermediate orientation, 
with the bar again along the x axis and a galaxy inclination of 77°, 
i.e. similar to that of M31, so as to allow a first comparison. In 
the innermost region the classical bulge dominates and gives near- 
elliptical isodensities. The extent and shape of this region depends 
on the mass and scale-length of the classical bulge compared to the 
disc mass and scale-length. Somewhat further out the boxy bulge 
dominates. This is a part of the bar, which when seen face-on is 
not distinguishable from the bar (see e.g. bottom panels of Fig. 0. 
When the galaxy is seen edge-on and the bar side-on (i.e. with the 
line of sight along its minor axis) the boxy bulge is seen to swell 
out of the equatorial plane (upper panels of Fig.0. When seen at 
an angle similar to that of M3 1 , this gives a more or less (depending 
on the model) boxy outline. For some models this is similar to the 
outline of the isophotes in the boxy bulge region in M31, as will be 
seen in Sect. l3.ll The fact that it is the boxy bulge that dominates 
the dynamics and sets the isophote characteristics in this region led 
us to term this isophotal region as boxy region. 

In regions further out than the boxy region the isophotal shape 
varies significantly from one model to another. For model MD the 
isophotes show no conspicuous feature and are rather similar to 
those in the boxy region. The isophotes in model MH show a very 
characteristic pinching toward the center from both sides of the bar 
minor axis (i.e. from above and below the galaxy major axis). This 
is similar in shape to the pinching that creates the peanut or the 'X' 



shape in the edge-on views and is in fact due to it. The part of the 
isophote that is near the bar major axis shows two strong protuber- 
ances on either side of the centre. The two remaining models (MHB 
and MDB) have in this region isophotes of similar shape, so they 
can be described together. Due to the high inclination and to the 
fact that the boxy bulge sticks well out of the equatorial plane, the 
boxy shape of the bulge will still be clearly visible near the minor 
axis. Near the major axis, however, there will be no contribution 
from the boxy bulge (because of its shorter extent) and one will see 
directly the outer part of the bar, which is vertically thin, i.e. does 
not extend much outside the equatorial plane. This contributes an 
extension, or elongation, of the isodensities towards the direction 
of the galaxy major axis and we will refer to the region where this 
is seen as the flat bar region. 



3 COMPARISON WITH THE NIR PHOTOMETRY OF 
M31 

3.1 Comparison of observed isophotes with model 
isodensities 

Fig.|2|displays the observed J image of M31, presented and anal- 
ysed in detail in Paper II. The isophotes show different characteris- 
tics in different regions of the galaxy. Starting from the center out- 
wards, we first have a center-most region with elliptical isophotes, 
presumably due to the classical bulge. We then encounter a region 
where the isophotes are boxy-like (see second and third isophotes 
from the centre in Fig.|2j- A similar region was seen in our mod- 
els in the previous section and we called it the boxy region. Note, 
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Figure 2. Grey scale representation and isophotes for the J image of M3 1 . 
North is at the top and East to the left. 

however, that contrary to the model, these boxy isophotes are skew 
with respect to the galaxy major axis, a feature that we will explain 
below. Beyond this region there is an abrupt change of the isopho- 
tal shape (see e.g. the fifth isophote from the centre in Fig.[2j. Here 
two of the sides of the isophotes (the ones near the galaxy minor 
axis) are quite straight, similar to the shape of the isophotes in the 
boxy region. However, near the major axis, which is at 38°, the 
isophotes have a clear elongation pointing to a direction somewhat, 
but not much, offset from the direction of the galaxy major axis. 
This is similar to what is seen in two of our fiducial models, namely 
the two models with a classical bulge. Thus, the shape of the M31 
isophotes constitutes a strong argument that M31 is barred. Let us 
now examine whether they will allow us to go further and to set 
constraints on the bar parameters. 

An important feature to note here is that the isophotes of M3 1 
do not show any pinching towards the center and this allows us to 
eliminate very strong bars, such as those in model MH, as likely 
candidates for reproducing the isophotal structure of M31. Indeed, 
the horizontal and vertical structures of a bar are strongly cou- 
pled. As shown by the simulations of A05 (see also AM02), the 
strongest bars have clear peanut (or even 'X') shapes when viewed 
edge-on and side-on, while the less strong bars have a boxy side-on 
shape. This is in good agreement with the observational results of 
iLtitticke et al. I l2000t) . who showed that stronger bars have more 
prominent peanut shapes. Peanuts, however, and even more so 'X' 
shapes, have a narrowing of the isophotes/isodensities towards the 
center if there is no bulge component and this is preserved also 
for viewing angles not far from edge-on. Thus, strong bars viewed 
at such angles would always show this characteristic squeezing or 
pinching of the isophotes in the central parts and the fact that this 
is absent in M31 argues that its bar can not be very strong. The 
isodensities of the MD model also do not reproduce well the ob- 
served isophotes, since, in the flat bar region, they miss the char- 
acteristic elongation in a direction near to the galaxy major axis. 
Thus, MD-type bars are also inconsistent with the M3 1 bar. In fact, 
only the two models with a classical bulge, MHB and MDB, give 
a good representation of the isophotes in all the relevant region. 
Even those, however, do not reproduce the slight skewness of the 
isophotes (i.e. their asymmetry with respect to the galaxy major 
axis) in the orientation shown in the middle panels of Fig.Q This, 
nevertheless, can be achieved if we consider a different viewing 
angle. 

Three angles are necessary in order to determine the orienta- 
tion of the M3 1 disc and bar, namely the inclination, the angle be- 



tween the bar major axis and the galaxy major axis and the position 
angle of the galaxy. This third angle does not contain any physical 
information, and will not affect the comparison of isodensity and 
isophotal shapes. The constraints brought by our simulations to the 
inclination angle are not tight. For one thing, there is no way of 
distinguishing between front and back, since A-body models do 
not contain dust. Moreover, the axial ratio of the inner isophotes is 
a function of three things : the inclination angle of the galaxy, the 
bar axial ratio in the galaxy equatorial plane (since a fatter bar will 
need to be viewed nearer to edge-on than a thinner one to give the 
same projected axial ratio) and the vertical extent of the bar ma- 
terial. Thus, the observed axial ratio in the inner regions can not 
determine uniquely the inclination angle. 

In Fig. [3] we vary, for model MHB, the angle between the 
bar and the galaxy major axis, to check its effect on the isoden- 
sity shapes. When this angle is zero, i.e. when the bar lies on the 
galaxy major axis, all isodensities are symmetric with respect to 
this common axis, as expected. This symmetry is broken as soon as 
this angle becomes nonzero. Then the boxy isophotes in the inner 
parts become skew, while the elongation of the isodensities some- 
what further out is offset from the galaxy major axis. Both these 
features are seen in the corresponding isophotes of M31, thus ar- 
guing that the bar in this galaxy is at an angle with the line of 
nodes. The direction of the skewness of the boxy isophotes and 
the direction at which the tips of the isophotes in the flat bar re- 
gion are pointing gives an indication about the position angle of the 
bar. Both directions are pointing further to the SE of M31's north- 
ern semi-major axis and this tells us that the position angle of the 
bar major axis is larger than the position angle of the galaxy. The 
amount of skewness and the corresponding isophotal shapes give us 
some constraints on the difference of these angles, as can be seen in 
Fig.|3| The best fits are for 20° and 30° (measured on the plane of 
the galaxy). Indeed, for 0° and 10° the asymmetries and skewness 
are insufficient, while angles of 40° and larger again do not give 
a good representation. It should be noted, however, that the exact 
values of this range may depend on the model. 

Fig. |4] uses appropriate orientations for all four models, to 
allow a full comparison with the M31 isophotes. It shows that, 
even with an optimum angle between the bar and the galaxy ma- 
jor axis, models MD and MH do not fare well. On the other hand, 
models MDB and MHB reproduce well the characteristic of the 
M3 1 isophotes, including the observed skewness. Thus, comparing 
isophotal and isodensity shapes has given us precious information 
on the existence and the properties of the M31 bar. 

3.2 Radial luminosity and density profiles 

In this section we compare radial luminosity profiles made along 
cuts on M31 as projected on the sky and radial density profiles ob- 
tained from the models in a similar way. 

Using the J image of Paper II, we obtained the radial lumi- 
nosity profile on a cut along the galaxy major axis, as well as on 
cuts parallel to it but offset by multiples of 50" either to the SE or 
to the NW. The results are given in Figs.|5]and|6| for the SE and the 
NW, respectively. In all profiles we note interesting asymmetries 
between the two sides. First, the position of the maximum drifts 
with respect to the center with increasing offset. For offsets to the 
NW, the maximum drifts towards the SW, while for offsets to the 
SE it drifts towards the NE. Note also that on the SE cuts there is 
a hump on the profiles to the NE of the center. In fact, one can see 
there a plateau ending at about 1000" from the nucleus, followed 
by a relatively steep drop starting at 1000" and ending roughly at 
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Figure 3. Isodensities for model MHB and various angles between the bar and the galaxy major axes. From left to right and top to bottom this angle varies 
from to 70° in steps of 10°. The inclination angle is in all cases 77°. 




Figure 4. The four fiducial models viewed in a M31-like orientation, i.e. an inclination of 77° and a position angle of 38°. Measured on the plane of the 
galaxy, the angle between the galaxy major axis and the bar major is 10° for the two models in the left panels and 20° for the two models in the right panels. 
The isodensities are chosen at levels showing best the features discussed in Section lTTl 



1300". The intriguing thing here is that this plateau occurs only on 
the one side of the centre. Such an asymmetry is also seen on the 
cuts parallel to the galaxy major axis and offset to the NW of the 
galaxy (Fig.[3J- In this case, however, the bump is on the SW side 
of the nucleus and is less prominent. 

The existence of the humps and the asymmetry of their lo- 
cation on the cuts provide useful clues about the existence and 
the properties of the bar in M 31. Indeed, the observations of 
lElmegreen & Elmegreenl 1 19851) show such humps in similar cuts 
in strongly barred early type barred galaxies. In that paper, the cuts 
were made along the bar major axis and show that the humps are 
symmetrically located with respect to the nucleus. Such humps are 
also seen for model MH in AM02 (see their Fig. 5) and they are 
also symmetrically located on either side of the center. On the other 
hand, no such humps are seen in most late type barred galaxies 
lElmegreen & Elmegreerll985h . where the decrease of the density 
is exponential-like. A similar decrease is found for model MD in 



AM02 (their Fig. 5). This could be a clue that the M31 bar is more 
likely to be MH-type than MD-type. 

Let us now follow the clues given by the asymmetric location 
of the bumps. Humps alone are not unambiguous evidence of the 
existence of a bar, since they can also be made by a lens-like struc- 
ture with abrupt edges, even if the lens is circular or near-circular. 
Such a structure, however, could not be responsible for the M31 
humps, since any axisymmetric, or near-axisymmetric, structure 
would place the humps symmetrically with respect to the nucleus. 
Thus, the existence and the asymmetry of the humps, when taken 
together, is an argument for the existence of a bar in M3 1 . If the ma- 
jor axis of this bar coincided with the galaxy major axis, the humps 
would again be located symmetrically with respect to the nucleus. 
Thus, the asymmetry of the bump locations argues that they are due 
to a strongly non-axisymmetric feature, i.e. a bar, whose major axis 
is at an angle with respect to the galaxy major axis. This result is 
in good agreement with was found in the last section by comparing 
isophotes to isodensities. 
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Figure 7. The upper panels show the luminous material of our four fiducial simulations, viewed at an inclination angle of 77°. For the two models on the left, 
the angle between the bar major axis and the galaxy major axis, measured on the plane of the galaxy, is 10°. For the two models on the right it is 20°. The 
horizontal stripes delineate the regions from which we calculate the density. The lower panels show the logarithm of the density as a function of position along 
these stripes. The solid line corresponds to the stripe outlined with the solid horizontal lines in the upper panels, and the dashed one to the stripe outlined with 
the dashed horizontal lines. The vertical lines give the lengths of the bar and of the boxy bulge, as found from Fig.Qif no projection is taken into account. 
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Figure 5. Luminosity of M31 along cuts parallel to the galaxy major axis 
and offset from it towards the SE by 0", 50", 100", 150", 200", 250" 
and 300". The distance along the cut is measured on the abscissa in arc 
seconds. The ordinate is labeled in mag / arcsec 2 , to within an arbitrary 
constant. Note the characteristic asymmetry of the luminosity profile with 
respect to the center. 



The location of the bumps can also give indications on the 
angle between the galaxy and the bar major axes, as well as on 
the length of the bar. We show this for our four fiducial models in 
Fig,|7J For models MH and MD we rotated the disc so that the left 
(eastern for comparison with observations) part of bar major axis 
is at an angle of 10° to the left (eastern) part of the galaxy major 
axis and then projected it to an inclination angle of 77°. We did the 
same thing with models MDB and MHB, except that now we chose 
the angle between the bar and the galaxy major axis to be 20°. We 
used 10° and 20° , respectively, because isophote to isodensity com- 
parisons (Sect. [33 show that these angles best reproduce the M31 
observations. Other values of these angles will be discussed later in 
this section. The results are given in the upper panels of Fig. We 
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Figure 6. As in Fig.|5] but for cuts offset towards the NW. 



then make cuts parallel to the galaxy major axis with offsets above 
and below. The density along these cuts can be seen in the lower 
panels of Fig.0 We find that the cuts for model MD are symmetric 
with respect to the center, the maxima in the two stripes coincide 
and there are no humps on either of the cuts. This argues strongly 
that an MD type model is a very unlikely candidate for the M31 
bar. The cuts for model MH show an asymmetry, the maxima in the 
two stripes are offset and there are humps asymmetrically located, 
as in M31. The amplitude of these humps, however, is too large 
and thus this model is not a good candidate for M3 1 . On the other 
hand, models MHB and MDB fare much better. Their profiles are 
asymmetric and in the same way as M3 1 . In particular, the position 
of the maximum shifts from the center in the same direction as for 
M31, i.e. towards the left (eastern) part for cuts to the South and 
towards the right (western) part for cuts to the North. Furthermore, 
we note the existence of a single hump on each profile, located in 
the direction towards which the maximum of the profile is shifted, 
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again in agreement with what is observed in M3 1 . This argues that, 
as in the MHB and MDB models in Fig.0 the eastern semi-major 
axis of the bar in M31 is to the South of the eastern galaxy semi- 
major axis. If it were to the North, then the asymmetries would be 
reversed. Thus, the sense of the asymmetries on the radial luminos- 
ity profiles give us constraints on the way the bar is oriented with 
respect to the galaxy major axis. 

Such cuts give us information not only on the sign of the angle, 
but also constraints as to its value. We repeated the above analysis 
considering different angles between the bar and the galaxy major 
axis in steps of 10° (results not shown here). For model MD we 
find that no angle gives a reasonable fit. For model MH, the best fit 
is for an angle of 10°, but, as can be seen in Fig.Q even that is not 
good. For model MDB we find that the best fit comes for an angle 
of 20°, but 30° is still good. Finally, for model MHB both 10° and 
20° are good. Thus, we can conclude that the most probable angle 
is about 20°. This is in good agreement with the range of values 
we found in the previous section. Adopting the M31 disc position 
angle of 38°, this means that the position angle of the M31 bar is 
of the order of 45°. 

Further information can be obtained on the length of the bar. 
For this we plotted in Fig. Q vertical lines at the radii at which we 
estimated the length of the bar and of the peanut component from 
the face-on and end-on views, respectively (Sect.|3- Note that the 
bar length is beyond the end of the flat part of the hump, nearer 
to the end of the abrupt drop after the plateaux. Using this, it is 
possible to make a rough estimate of the length of the bar in M31. 
As seen in Fig. [5] the flat part of the bump there extends till 1000" 
from the center, while the abrupt drop after the flat part extends 
from 1000" to 1300". Thus, we can estimate the length of the bar 
in M31 to be of the order of 1300". Of course this estimate has 
a considerable uncertainty, since we can not determine the end of 
the steep drop in the profile precisely either in the model, or in the 
data. There is a further uncertainty from the fact that the bar semi- 
major axis is at an angle to the galaxy semi-major axis, which we 
have only roughly constrained. Changes of the order of 2%, 5% or 
14% are to be taken into account if this angle is 10°, 20° or 30°, 
respectively. 

For model MH we plotted two vertical lines for the bar length, 
corresponding to the two estimates of the bar length given in Fig.Q 
The shorter of the two is obtained if we consider that the bar ends 
where the ring starts and the longer one if we assume that the bar 
continues within the ring reaching its outer end. From Fig. it is 
clear that the end of the sharp drop after the end of the plateau 
coincides with the second definition, giving the longer barlength. 
This, however, does not mean that it is this definition that is correct, 
since clearly the ring participates in the formation of the hump. 
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Figure 8. Rotational velocities of M31 from Rubin & Ford (1970). Note the 
deep minimum around 10'. 




afl'S s Mrs 



Figure 9. Observed PVDs from the HI in M31 from cuts parallel to and 
near the galaxy major axis. Reproduced from Brinks & Shane (1984). 



4 KINEMATICS 

The kinematics of M3 1 have been well observed, both in the HI and 
the Ha. Our A-body models, however, do not contain gas, so it is 
not possible to make any direct comparisons to the corresponding 
observations. Nevertheless, the observed gaseous kinematics give 
strong arguments for the existence of a bar in M31. For this rea- 
son, we will discuss here first the relevant observational results and 
then present simulated P VDs of the gas compo nent in a strongly 
barred galaxy model from I AthanassoulJ <1992d) . Although this is 
an idealised model (with a Ferrers H 8771 ellipsoidal bar) and not a 
simulation result as the four fiducial models we discuss here, it will, 



nevertheless, allow us to discuss some generic features of PVDs in 
barred galaxies. 

The first clear kinematical clue that M31 is barred was con- 
tained in the data of iRubin & Fordl <1970h but, at the time, was not 
interpreted as such. The velocity measurements from that paper, 
read off their Fig. 9, are reproduced in Fig. [8] They show a very 
deep minimum near 10', which is not observed normally in spiral 
galaxy rotation cu rves. This min imum, described as 'the annoying 
deep minimum 'bv lRubinlll994l) . provoked discussion, often going 
as far as to cast doubts as to the quality of the data. 

Brinks & Shane (1984) observed M31 in HI and gave a num- 




Figure 10. Response of the gas to a bar. The major axis of the bar is at 45° 
to the horizontal and is rotating clockwise. It has a semi-major axis length 
of 5 and a semi-minor axis of 2 length units. Lighter shades denote higher 
density regions. 

ber of PVDs on cuts parallel, or perpendicular to the major axis. 
The cuts parallel to and near the major axis (some of which are 
reproduced in Fig. [5} show very interesting structure. In particular, 
they show two branches, separated by a lo w density region. In order 
to explain and reproduce these features, iBrinks & Burton! Il984l) 
proposed a model of a rotating disc which flares and warps in its 
outer regions. According to their model, the inner branch would be 
due to gas near the galaxy major axis, which would imply that the 
rotation curve for that galaxy would be very sharply rising. It would 
then be the outermost gas which, projected onto the line of sight, 
is responsible for the lower branch (the one where the velocities 
rise slower with distance from the center). Their model reproduces 
a number of the main characteristics of the PVD diagrams (namely 
the existence of two branches through the center), but can not be 
qual ified as a good fit to the data. 

iRubinl J 1994) su perposed the iRubin & Ford Jl97d) velocity 
measurements on the IBrinks & Shand PVDs going through the 
center of the galaxy and found good agreem ent. She thus e stab- 
lished that there was nothing wrong with the IRubin & Ford data 
and stressed that the low velocities should not be identified as ro- 
tational velocities. She also proposed that the complex velocities 
could be due to extra-planar gas in the innermost parts, without, 
however, showing that any such model could inde ed reproduce the 
PVD characteristics. Thus, both lBrinks & Burton! and lkubinl relied 
on extra-planar gas in the central and/or in the outer parts of the 
galaxy, but did not invoke the possibility of a bar. It was only a few 
years later that it became clear that such structures arise naturally 
if the galaxy is barred. 

Structure such as seen in the above discussed PVDs has 
been witnessed in the observed emission line PVDs of edge-on 
galaxies with boxy or peanut bulges and was readily attributed 
to the bar jKuiiken & MerrifielJl995tlMerrifield & Kuiikerll999t 
iBureau & Freemadl 19991) . Indeed, gas responding to a strong bar 
potential does not fill the bar region uniformly, but, on the con- 
trary, forms strong c oncentrations and relatively empty regions 
lAthanassouklll992bh . This is illustrated in Fig. 1101 which shows 
the gas response to a strong bar model (model 047 of lAthanassoulal 
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ll992df ). We see that most of the bar region is relatively empty 
of gas except for the region surrounding the center and two nar- 
row stripes along the leading ed ges of the bar. These are in fact 
the loci of shocks in the gas flow l Athanassoula 1992b). Secondary 
concentrations can be seen sur rounding the bar region and at the 
ends of the bar. As shown by lAthanassoula & Buread 1 19991) . it 
is the existence of these high density and low density regions, to- 
gether with the corresponding velocities, that creates the charac- 
teristic structures of the PVD s in barred galaxies. To show this, 
lAthanassoula & Buread 1 19991) obtained model PVDs by project- 
ing edg e-on some of the two-d imensional hydrodynamic simula- 
tions of I Athanassould 1 1992b]) . We repeat this for the model in 
Fig. llOl and display the result in Fig. II II One can immediately see 
the two branches that stand out in the M31 PVDs, separated by a 
low density region. The inner branch is due to the gas near the cen- 
ter, where the flow lines are elongated roughly along the bar minor 
axis. The outer and lower of the two branches is due to material in 
the outer regions. The two branches are separated by low density 
regions on the PVDs, because the two regions are separated by low 
density regions in the face-on view of the bar. 

Further comparison between Figs.l9land ll H allows us to check 
whether the values for the bar position angle and length, as obtained 
from the photometry, are consistent with the observed PVDs. In the 
previous sections we found the best fit between models and obser- 
vations when the angle of the bar major axis is around 20° from 
that of the galaxy major axis. This is near the viewing angle of the 
upper middle panel in Fig. II II We note that in this case the contri- 
bution from the central part, i.e. the inner branch of the PVD, is a 
nearly straight segment passing through the center of the PVD, in 
agreement with the observations. Also, this segment rises to veloc- 
ities higher than those of the outer parts, again in agreement with 
the observations 5 . Note that this latter check excludes a large value 
of the angle between the bar and galaxy major axes, as e.g. in the 
lower middle panel of Fig. II II A further consistency check can be 
made with the bar length. From Fig. II H it can be seen that the max- 
imum extent of the low density region between the two branches of 
the PVD is somewhat smaller than the length of the bar semi-major 
axis. Applying this to Fig.|9|we see that the barlength of 22', which 
we estimated in Sect. l3.2l based on the photometry, agrees with the 
observed PVDs. 

Unfortunately, the analysis of the kinematics can not be 
pushed any further than the discussion of generic bar features for 
many reasons. The bar model used in this section is an idealised 
bar model and not one of the four fiducial models discussed in the 
previous sections. Moreover, the PVDs of Fig. fTTI were obtained by 
integrating along the line of sight a perfectly edge-on razor thin gas 
layer, whereas the inclination angle of M31 is 77°, its gas layer is 
not razor-thin and this layer could also flare and warp in the out- 
ermost parts. Although this is far outside of the established bar re- 
gion, it could still influence the PVDs if gas was folded onto the 
line of sight. This would give furthe r branches on the PVD s. In- 
deed, the M31 PVDs (see Fig. 8 of IBrinks & Shanelll984l) show 
suggestions of many more features which, due to the lack of gas in 
extended areas of M3 1 , can not be followed clearly. For example, 
in many of the cuts one has the distinct impression that there are 



5 As explained in Athanassoula & Bureau] ll999l) . this is due to the fact 
that in the central regions the gas flow lines follow roughly the orientation 
of the X2 family. The largest velocities occur of course near pericenter and 
would, in this region, be observed when the bar is viewed near side-on, as 
is the case. 
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Figure 11. Model PVDs for the gaseous component in the strongly barred galaxy used in Fig. 1101 Results are presented for six values of the angle between the 
galaxy and the bar major axis (given in the upper right corner of each panel). 



three branches. To model all these features fully is an interesting, 
but daunting task, further complicated by the lack of gas in some 
crucial regions, and well beyond the scope of this paper. 



The stellar velocity field of our four models is given in Fig. l 1 21 
Unfortunately, no corresponding observational data is available, so 
these figures can only be considered as predictions. For the viewing 
angles we adopted the values found to be optimum in the previous 
sections, i.e. an angle between the bar and galaxy major axes of 10° 
for models MD and MH and to 20° for models MDB and MHB. 
The inclination angle was taken to be 77° in all cases. When calcu- 
lating these velocity fields we included both the disc and, whenever 
relevant, the classical bulge material. This makes the effect of the 
bar less strong, but will allow comparisons with future observa- 
tions. 



The velocity field of model MD shows no clear-cut perturba- 
tions which could serve as clues to the existence of the bar. On the 
other hand, the isovelocities of model MH have sharp deviations 
due to the bar, which give t hem the form of a lette r 'Z', as seen 
in many barred galaxies (e.g. lPeterson et al\ ll978lV ). Such devia- 
tions, but less clear cut, are seen also for models MDB and MHB. 
A telltale sign for the existence of a bar or of an oval is the fact that 
the kinemat ical and photometric minor axes do not align iBosmal 
Il978tll98lh . This would have been easily spotted for models MH, 
MDB and MHB, as well as in a very careful analysis of the velocity 
field of model MD. 
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Figure 12. Stellar velocity field of the four models discussed in this paper, 
namely, from top to bottom, models MD. MH. MDB and MHB. The isove- 
locities are given by thick lines and the kinematic major axis with a dashed 
line. We also overlay the isodensities with thin contours. 
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5 DISCUSSION 

5.1 Comparison with previous work 

Look ing at the isophotes of M3 1 <LindbladlT95tl IdeVaucouleur i 
Il958l etc) one imme diately sees that the part considered as the 
bar in previous work <Lindbladll956tlStarkll977l:IStark & Binnevl 
1 1994) is in fact the boxy/peanut bulge and in the last two papers 
was often termed as such. Thus, the 'bar' they propose should be 
shorter than ours, and also its angle with the northern semi-major 
axis of M3 1 should be larger. Indeed, this is the case. In our model 
the bar semi-major axis length is about 22', the angle between the 
galaxy and the bar semi-major axes is about 20° and the axial ra- 
tio in the equatorial plane larger than 2. In the previous models the 
'bar' semi-major axis length is about 15' — 16', the angle between 
the major axes of the galaxy and the 'bar' is about 70° and the 
axial ratio in the equatorial plane smaller than 2. These latter num- 
bers are in agreement with those of the boxy bulge of our model. 
A more precise comparison is not possible since the 2D shape of 
our boxy/peanu t bulge is rect ang ular-like, while that ass umed by 
Lindblad il956l) . IStarkl J 19771) and lStark & Binnevl 1 1994) is an el- 
lipsoid. 

5.2 Spiral Structure 

The most prominent feature in M31 is a ring around 50'. More 
precisely, this is a pseudo-ring, since it has an opening in the SW 
and its distance from the center is a function of the azimuth. On top 
of the bar we are advocating here, M3 1 has two companions (M32 
and NGC 205), both of which can drive spirals and rings. 

M32 (NGC 221) is a small companion in the SE quadrant of 
M31. Baade believed it to be below the plane of M31 and its sys- 
temic velocity is compatible with a near-circular orbit which is ret- 
rograde with respect to the sense of rotation of M3 1 . Such an orbit 
excites an m = 1 inner Lindblad resonance in the target disc and 
thus drives a 1 -armed leading retr ograde spiral, as was shown ana- 
lytically by Athanassoula ( 1978 a]) and confirmed by iV-body sim- 
ulations by rThomasson et al I Jl98gf) . iBvrd & Howard i 19921) and 
IVozikis & CaranicolasHl994) . It can also be understood intuitively 
since a retrograde driving can not excite any 2-armed spiral 6 . So, 
if the motion of M32 is indeed retrograde with respect to the sense 
of rotation in M31, one would expect to see a very tightly wound 
1-armed leading spiral, comparable to a pseudo-ring, in a region 
compatible with the position and velocity of M32. Such a spiral 
component can indeed be found from the distrib ution of HII re - 
gions and OB assoc iatio ns, as has been shown by |Kalnajs| l 19751) . 
ISimien et al. I 1 19781) and|Considere & Athanassoulal 119821) . This, 
together with other comparisons with the morphology and kine- 
matics of M31, show that this is a plausible exp lanation to the 
M31 ring-like structure at a radius of ro ug My 50' <Kalnais|[l975l 
lAthanassoul Jl978bllslmien et al. Il97^> . 

Braun (1991) describes the spatial distribution of the HI gas in 
M3 1 by a global two-armed trailing spiral and obtains a good fit if 
he allows variations of the orientation of the galaxy plane both with 
radius and azimuth (i.e. local line of nodes). He notes that these 
mid-plane departures lack bisymmetry and also make the existence 

6 M32 has a mass of only a few percent of that of M3 1 , so that, in order to 
drive a spiral structure whose amplitude is compatible with that observed, 
it needs to rely on a resonance. A simple back-of-the-envelope calculation 
shows that a companion on a retrograde orbit can not induce an m=2 reso- 
nance, only an m=l one. 



of a massive dark halo unnecessary. He mentions two possibilities 
for the driving of the two-armed spiral, name ly the bar-like triaxial 
bulge <Starkll977llHodge & Kemiicuttll982l ; etc) and the nearby 
companion M32. He favours the second alternative due to the sub- 
stantial departures from a planar gas distribution. Since he does 
not give a dynamical model providing information on whether the 
mass and velocity of M3 1 are compatible with the observed spiral 
feat ures, it is not possible to assess this further. 

iGordon et al. I l2006t) also propose a model relying on M32, 
but without invoking any resonance. This, however, produces only 
a very large number of arms (m of the order of 10, see the central 
right panel of their Fig. 3) and no ring-like structure, i.e. a structure 
very unlike what is observed. Furthe rmore, if they use a mass of 
M32 in agreement with observations (Mateo 1998J), they find that 
'very little effect is seen in the M31 disc', in agreement with what 
could be foreseen from the simple explanations given above. They, 
thus, artificially increase the mass of M32 by a factor of 5 with 
respect to Mateo's estimate, suggesting that M32 could be more 
massive than originally thought. 

How does our proposed bar fare with the observed pseudo- 
ring? Taking for the semi-major axis length of the bar a value of 
at — 22', as proposed in Sect. 13.21 we can set constraints on 
the location of corotatio n, which should necessarily be outside the 
bar jContopoulosll98rJ) . Hydrodynamic simulations lAfhanassoulal 
Il992bl) have given a range of (1.2 ± 0.2)at for the co rotation ra- 
dius and this agrees well with observational predictions lElmegreenl 
ll99dlGerssen. Kuiiken & Merrifieldl2003l etc). This places coro- 
tation in the range 26' ± 4'. Assuming a flat rotation we then will 
have the outer Lindblad resonance of this pattern at 45' ± 7'. It 
should be noted that this error bar is based only on the range found 
from the hydrodynamical simulations, while, if we we took into 
account uncertainties on the end of the bar and on the form of the 
rotation curve, the error bar would be considerably larger. Thus, if 
the rot ation curve in the relev ant region is slightly rising (as shown 
e.g. in ICarignan et al. II2006T) . the location of the outer Lindblad 
resonance would be somewhat further out. For example, assuming 
that the velocity increases roughly as r° 1 places outer Lindblad at 
49' ± 8'. Barred galaxies are known to have outer rings and both 
analytical work (e.g. |Sjjiwarz|1198 1 ) and statistical analysis of the 
ring radii lAthanassoula et al. II 19821 : iButall 1 995f) argue that these 
are placed at the outer Lindblad resonance. Given the uncertainties, 
we can say that the location of the ring in M31 is consistent with 
its being an outer ring due to the bar. 

The two alternative explanations leading to a reasonable re- 
production of the M3 1 ring are thus both linked to resonances. In 
the first one, the pseudo-ring is a very tightly wound one-armed 
spiral formed by a resonance with M32 rotating in a retrograde or- 
bit around M3 1 . In the second one, it would be formed at the outer 
Lindblad resonance of the bar. Is there any way of discriminating 
between the two alternatives? An argument in favour of the first 
explanation is that the observed ring is a pseudo-ring and its de- 
viations from a ring are more of the m = 1 than of the m = 2 
type. This, however, is not a strong argument, since these devia- 
tions could be due to the effect of M32 on a pre-existing closed 
outer ring. An argument in favour of the second explanation is that 
it is statistically more likely, and would make M3 1 a more ordinary 
galaxy, its pseudo-ring in good agreement with a structure com- 
monly found in a large fraction of barred galaxies. Of course this is 
not on its own a sufficient reason to eliminate the first explanation, 
particularly since, if M31 was at a larger distance and observed 
with lower resolution, the tightly wound one-armed spiral would 
be confused with a standard outer ring so that several such struc- 
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tures could be actually observed and mis-interpreted. A more in 
depth modeling of M31 would be necessary to distinguish between 
the two possible explanations, or, better still, to produce a model in 
which both the effects of the bar and of the companions are taken 
into account. At this stage, we can only say that the second expla- 
nation is statistically more likely, but we still need to model the 
effect of M32 on the ring to see whether it is compatible with the 
observations. 

5.3 Relative lengths of the bar and of the boxy bulge 

Because of its viewing angle (near to but not quite edge- 
on), M31 brings crucial input to the issue of the rela- 
tive bar and boxy bulge lengths. Both orbital structure stud- 
ies and iV-body simulations argue that the extent of the 
boxy (or equivalently the peanut, or 'X') structure is shorter 
than that of the bar jSkokos, Patsis & Athanassoulal 120021 : 
IPatsis. Skokos & AthanassoiiiaM^OOZ A05). Observational confir- 
mation, however, is non-trivial, since the length of the bar can only 
be measured in near face-on galaxies and the length of the boxy 
bulge only in edge-on galaxies. The main input thus comes from 
edge-on galaxies on which the extent of the boxy bulge can be mea- 
sured directly, while the length of the bar is inferred from the length 
of th e plateau of the light p rofile on a cut along the equatorial plane 
(e.g. iLuttickeef aTlbOOOt AM02; iBureau & AthanassoulJl2005t 
A05). This measurement, however, is not very accurate, since the 
end of the plateau tapers off into the disc. It, nevertheless, provided 
clear evidence that the bar extends further out than the boxy bulge. 
M3 1 allows for a different way of comparing the two extents. Mea- 
suring along the major axis, one can delineate the region where the 
isophotes are boxy-like and the region where they have elongations 
pointing along the bar major axis, i.e. the regions named in Sect.|2| 
the boxy and the flat bar regions. The extent of these two regions 
provides us with the length of the boxy bulge and the length of the 
bar, respectively. This is well illustrated in Fig.Q Here the length 
of the bar is obtained from the face-on views (lower panels) and 
shown with the full vertical solid lines, which are continued on the 
two other panels to show where the end of the bar lies in the edge- 
on and M31-like viewing angles. The length of the boxy bulge is 
obtained from the edge-on view (upper panels) and plotted with the 
dashed vertical lines, which are continued on the two other panels 
to show where the end of the boxy bulge lies in the face-on and 
M31-like viewing angles. Looking at the middle panel we see that 
the extent of the regions with the boxy isophotes and the extent of 
the region with the elongation towards the bar major axis indeed 
provide us with the boxy bulge length and the bar length, respec- 
tively. 

Unfortunately, these lengths are not easy to measure directly 
from the isophotes. A better determination can be obtained from 
isophotal fits and calculation of the cn and &4 coefficients along 
them. Such fits will be made and discussed in Paper II. The mea- 
surements we have here, although not as accurate, are still sufficient 
for inferring the type of periodic orbit families which are involved 
in the building of M31's box. Using for bar semi-major axis length 
th e value of 22' and for th e boxy bulge extent the value calculated 
bv lStark& Binnevl < 19941) . i.e. 15' gives a ratio of 1.4. Although 
this number has a large uncertainty, as discussed in the previous 
sections, it seems safe to say that it is less than 2. This limit is 
important since it gives us indications that the main families build- 
ing the peanut are those linked to the higher order vertical reso- 
nan ces (e.g. families xlv3, or xlv4) and excludes the xlvl fam- 
ily IPatsis. Skokos & Athanassoulall2002T) . It would be important 



to repeat this calculation for a number of galaxies seen at roughly 
the same orientation as M31 and showing similar features, to see 
whether the higher order families are always the most probable 
building blocks for boxy bulges. 

5.4 Miscellanea 

The comparisons of simulations to observations argue that M31 has 
two types of bulges, a boxy/peanut bulge and a classical bulge. In 
that it is not unique. In fact, judging from the fact that we see in 
many far-from-edge-on early type barred galaxies both a classical 
bulge and a strong bar (indicating the presence of a boxy/peanut 
bulge), a large fraction of early type barred galaxies should have 
the two types of bulges. In late types, there should be many cases 
with both a discy bulge (A05) and a boxy/peanut bulge, but the 
fraction will be more difficult to assess. Finally, in many galaxies 
all three types of bulges could coexist. 

Our two 'best' models presented here are not meant to be an 
exact fit to M31. We presented four fiducial models to argue that 
M31 is barred and to set constraints on its bar properties. We have 
thus outlined a class of models which would give reasonable fits. 
Surely, by looking through a large number of models within that 
class it would have been possible to find some that gave better fits. 
This task, however, does not seem warranted for two reasons in 
particular. First, the free parameter space is very large. Indeed, lit- 
tle is known on the mass distribution of the halo (shape and radial 
profile) and even less about its velocity distribution. Even for the 
visible material, there is a lot of freedom in selecting Q(r). Fur- 
thermore, M31 has two companions, which will induce asymme- 
tries which any thorough modeling should also take into account. 
Thu s, a full mode l for M 3 1 is well beyond the scope of this paper. 

using B, V, R and I photometry of the 
center-most parts of M3 1, revealed a small inner bar of semi-major 
axis length about 0/5. Its axial ratio is approximately 3 and it 
has dark lanes on its leading edges which are slightly curving in 
the trailing sense. This small bar has about the same colour and, 
therefore the same population as the bulge stars. Such inner bars 
are found in the innermost parts of rough ly 30% to 40% of bars 
jLaine et al. l2002UErwin & Sparkj200a . 

The local group contains three large spiral galaxies - namely 
our own Galaxy, M31 and M33 - and all three were initially con- 
sidered as unbarred, de Vaucouleurs, based on the so-called for- 
bidden gas velocities towards the galactic center, proposed first in 
1964 that our galaxy was barred and th is is now consid ered as a 
well established fact (for a review see |Pehnenll2002bl and ref- 
erenc es therein). The first arguments that M31 is barred came in 
1956 lLindbladll956l) . and we have presented in this paper a num- 
ber of strong arguments supporting this. Finally, images in the red, 
or NIR of M33 show that the spiral arms in that galaxy do not reach 
the center and that in between them there is indication of a small 
bar. It thus looks like the two, or perhaps even the three, largest spi- 
rals in our local group are barred. This is in good agreement with 
the results found in the NIR from larger samples of disc g alaxies 
lEskridge et a71l200Ct lGrosb0l et g/"l2004 iKnapenlli^l . show- 
ing that the vast majority of spirals is barred, although in some 
cases the bar is small. 

The isophotal shapes of M31 are very characteristic of boxy 
galaxies, providing straightforward evidence for the existence of 
the bar. As previously discussed, NGC 7582 and NGC 4442 have 
a similar orientation. In particular, the isophotes of the latter 
iBettoni & Gallettal 19941) have characteristics very similar to those 
of M31 (see their Fig. 2). As in M31, we can distinguish there the 



Unravelling the mystery of the M31 bar 13 



boxy and the flat bar region, as well as a similar skewness. Since 
these features are easily identified, it would be important to study 
a number of galaxies seen at roughly the same orientation as M3 1 
and showing similar characteristics, in order to get more informa- 
tion on the relative extents of the peanut and the bar and from that 
infer which orbital families constitute the building blocks for boxy 
bulges. 



6 SUMMARY 

In this paper we argue that M3 1 is a barred galaxy. Its bar is not 
as easily distinguished as it is in galaxies nearer to face-on, but it 
still leaves a number of strong clues, both in the photometry and 
in the kinematics. Our arguments are based on comparisons with 
TV-body simulations. Rather than using only one simulation, which 
would have allowed us to point out only the generic bar properties 
in M31, we use four models with a range of bar properties. In this 
way we can argue that there is a bar, as well as try to constrain 
the properties of this bar. Our four models are called MD, MH, 
MDB and MHB, respectively, for continuity with previous works. 
MD has the weakest of the four bars and MH the strongest. MDB 
and MHB have a classical bulge as well, contrary to the other two 
models which have only the boxy version. 

Comparison between the isodensity curves from our models, 
projected in the same orientations as M31, and the M31 isophotes 
argues very strongly for the existence of a bar. It also allows us to 
say that MH is a very unlikely candidate for the M3 1 bar, since this 
model has a characteristic pinching of the isophotes in the central 
parts, totally absent from the observed isophotes of M31. The MD 
model is also a rather unlikely candidate, because its isodensities do 
not show the required shapes. On the other hand, the isodensities of 
our MDB and MHB models reproduce well the observed isophotes. 
Both show two distinct regions with characteristic isophotal shapes. 
The inner of the two regions is dominated in the models by the boxy 
bulge and so we call it the boxy region. From the similarity of the 
isodensities and the isophotes in this region we can infer that in 
M31 also it is dominated by the boxy bulge. Outside this region ex- 
tends in the models the flat (vertically thin) bar, which gives to the 
isodensities characteristic elongations towards a direction near, but 
somewhat offset from the galaxy major axis. We call this region 
the flat bar region. Again these isophotal shapes are found in the 
M31 isophotes and we can thus attribute them to a bar. The M31 
isodensities show also a characteristic skewness, which can be well 
reproduced by our MDB and MHB models provided the angle be- 
tween the bar and the galaxy major axes is within a given range of 
values. 

We also compared radial luminosity profiles from M3 1 - made 
along slits parallel to the galaxy major axis and offset either to the 
SE, or to the NW - to radial density profiles obtained from the mod- 
els in the same way. Here again the bar has left a number of clues. 
The humps and asymmetries on the observed luminosity profiles 
argue clearly for the existence of a strongly non-axisymmetric com- 
ponent, such as a bar. Comparisons with the A^-body results allow 
us to constrain the models. Model MD can not reproduce observed 
luminosity profiles, while model MH reproduces them badly. It is 
again models MDB and MHB that give good fits. These compar- 
isons also give constraints on the value of the angle between the 
galaxy and the bar major axes and the range thus obtained agrees 
well with that obtained with the help of the isophotes. Furthermore, 
these comparisons allow to give a rough estimate of the length of 



the bar major axis, which we can estimate to be of the order of 
1300". 

The PVDs observed in the HI show also characteristic bar sig- 
natures. In particular there are two branches separated by an empty 
space. The inner branch is due to gas in the central region and the 
second one to gas outside the main bar region. The empty space is 
due to the lack of gas in the main bar region. These PVD diagrams 
also set constraints on the bar position angle and length. Both are 
in good agreement with the values found from the photometry. 

From all these comparisons we can conclude that there exists 
very strong evidence that M31 is barred. In our two best models 
the bar is neither too strong, nor too weak and its major axis has a 
position angle of roughly 45°, i.e. there is angle of about 20° be- 
tween the galaxy and the bar semimajor axes in the plane of the 
galaxy. These models have both a classical and a boxy bulge. The 
length of the latter is smaller than that of the bar and the most prob- 
able backbone for the M31 boxy bulge are periodic families from 
families of the xi-tree, bifurcating from the xi family at a higher 
vertical resonance, such as the the xlv3 or xlv4. A bar such as that 
of our models could create an outer ring, whose location would be 
in agreement with the pseudo-ring observed in M3 1 at roughly 50' . 
This gives a second, alternative, explanation for the pseudo-ring. 
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